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The internal symmetry of composite relativistic systems is discussed. It is demonstrated that Lorentz-Poincare 
symmetry implies the existence of internal moments associated with the Lorentz boost, which are Laplace- 
Runge-Lenz (LRL) vectors. The LRL symmetry is thus found to be the internal symmetry universally 
associated with the global Lorentz transformations, in much the same way as internal spatial rotations are 
associated with global spatial rotations. Two applications are included, for an interacting 2-body system 
and for an interaction-free many-body system of particles. The issue of localizability of the relativistic CM 
(^ . coordinate is also discussed. 
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p ,: I. INTRODUCTION 



Separation of the internal dynamics of composite relativistic systems from their global dynamics is one of the long- 
standing yet not fully (or satisfactorily) solved issues in fundamental theoretical physics. Lorentz-Poincare symmetry 
implies the constancy of the total linear momentum P^ and total angular momentum J^'^ , but not the uniqueness of 
the centre-of-mass (CM) coordinate. Thus, various definitions have been suggested for the relativistic CM, based on 
different requirements (see, e.g., Refs.[l|-yfor early publications on the subject, or Refs.Qandyfor recent publications 
^ ■■ with an extensive bibliography covering its history). 
OO In non-relativistic systems, the separation of the total dynamics into CM motion and internal dynamics is incorpo- 

T-H rated just in splitting the total energy and the total angular momentum into sums of CM-terms and internal terms. In 

OC ' relativistic systems we know, correspondingly, to isolate the internal energy as the total invariant mass of the system, 
" . , and to extract from J^'' its part which is responsible for internal rotations - the spatial internal angular momentum 
I/"} ■ tensor (see Eq. ([5]) below). 

^D , However, internal symmetry is more than just rotations. It is well known that in 2-body Newtonian Kepler- 

^^ ■ Coulomb systems there exists the Laplace-Runge-Lenz (LRL) vector as a constant of the motion. Knowledge of the 
'. . ' LRL vector amounts to having a full solution for the configuration of the system (details of orbit, etc.)*' ; in the 
corresponding quantum systems the LRL vector provides a very elegant means for obtaining the full quantum picture 
of the system (as in the case of the hydrogen atom)i. Although still regarded by many as corresponding to 'accidental' 
or 'hidden' symmetry, particular only to l/r potentials, it is known already for a long time that LRL vectors exist 
C^ ' in all rotationally symmetric systems^""—, including various relativistic models as well^i^^""— . As in classical Kepler- 
Coulomb systems, these LRL vectors may be constructed to be constant^^, and they generate, together with the 
internal angular momentum, 5*0(4) or 5*0(3, 1) symmetry according to the energetical state of the system. It follows, 
therefore, that the mere existence of the LRL vectors is a direct and essential consequence of the existence of internal 
rotational symmetry, and the 50(4) or 50(3, 1) symmetry thus generated is a natural extension of the rotational 
symmetry^"— i2^. 

Considering composite relativistic systems, since Lorentz-Poincare symmetry implies internal rotational symmetry, 
the foregoing argument then leads to expect the appearance of LRL vectors in all relativistic systems. In fact, not 
only do LRL vectors exist in relativistic systems (and therefore participate in generating the internal symmetry there) , 
their very origin lies within the relativistic domain : A number of years ago it was discovered by Dahl^^i^i that the 
Newtonian LRL vector emerges in the computation of the Lorentz boost in the post-Newtonian approximation of 
electromagnetic or gravitational 2-body systems. Although it is only the Newtonian LRL vector this an essentially 
relativistic result, because it appears in terms of order 1/c^, vanishing in the full non-relativistic limit when the 
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Lorentz boost becomes the Galilei boost. Dahl's results were recently re-established and extended to other systems, 
including fully relativistic non-interacting 2-body systems^ and arbitrary post- Newtonian centrally symmetric 2-body 
systems22i. 

Combining these two characteristics - the existence of LRL vectors for all rotationally symmetric systems together 
with their relativistic origin - it is then natural to expect that LRL vectors form an integral part of the internal 
relativistic symmetry. Also, since the spatial part of the relativistic CM coordinate is derived from the Lorentz boost, 
from which also the LRL vector is derived, it is natural to expect that the LRL vector plays a significant role in the 
determination of the relativistic CM. These subjects which so far were regarded as completely distinct - determination 
of the relativistic centre-of-mass, relativistic internal symmetry and LRL vectors - appear therefore to form part of 
one and the same story. 

So far, LRL vectors and the associated symmetry were considered, even in systems with relativistic features, as 
extensions of Newtonian systems rather than from the stand-point of full Lorentz- Poincare symmetry. The purpose of 
the present paper is thus to discuss in detail the internal symmetry in composite Lorentz- Poincare symmetric systems, 
in a manifestly covariant manner in Minkowski space-time, tracing the appearance of LRL vectors in these systems, all 
in association with the issue of determining the CM coordinate and its properties. The paper extends and completes 
recent publications on the subject^i^. 

The paper is composed in the following of three main parts. First, in Sees. HIl IIIII and IIVI are discussed in generic 
terms the appearance of the LRL vectors and their association with the CM coordinate. It is shown that the mere 
existence of Lorentz- Poincare symmetry implies the existence of LRL vectors, which together with the internal angular 
momentum generate the internal symmetry of the system, at least for 2-body systems. The spatial component of 
the CM coordinate contains an internal part which is proportional to the LRL vector of the system, and the type of 
the symmetry depends only on the energetical state of the system, being 50(4) or 50(3, 1) for bound or unbound 
systems, respectively. 

The methods developed in the first part are then applied in Sees. |V] and IVII to two particular systems. As an 
illustration of the application of these methods, a 2-body system with light-like antisymmetric scalar- vector interaction 
is considered in Sec. |Vl This is a relatively simple model, so chosen to keep the exposition simple and clear. Then, 
since the procedure developed in the first part allows the definition of many-body LRL vectors, it is applied in Sec. IVII 
to non-interacting many-body systems, thereby offering a first step into the generalization of the LRL symmetry also 
to interacting many-body systems. 

Finally, the implications of the foregoing results for the relativistic CM coordinate, in particular the issue of its 
localizability and its association with properties of the LRL vector of the system, are discussed in section I VIII Final 
discussion and concluding remarks are given in Sec. I Villi 

Notation. In the following we consider dynamics described in a Minkowski space-time {a;'^}, fi = 0,1,2,3 with 
metric tensor g^^ — diag(— 1, 1, 1, 1). Spatial (3D) vectors are denoted by an arrow. Dot product is A ■ B — A^B^. 
The unit fully anti-symmetric (Levi-Civita) pseudo-tensor is e^'^^f = —e^^xp = 1 for {jivXp) an even permutation 
of (0,1,2,3). It is also assumed throughout that c—\ unless specified otherwise. An orthogonality subscript (_L) 
implies the component of a 4- vector perpendicular to the total linear momentum P''. 

II. THE CM COORDINATE AND DECOMPOSITION OF THE TOTAL ANGULAR MOMENTUM 

We start by discussing the construction of the CM coordinate of general classical (non-quantum) composite rela- 
tivistic systems and the decomposition of the total angular momentum into CM-terms and internal terms. We consider 
systems of N point particles, with masses vria and moving on the trajectories x^ — 2:^(Ta), a = 1, ..., A'^, Tq being the 
proper time of the a-th particle. 

With constant total linear momentum, the space-time trajectory of the CM coordinate is expected to be a straight 
line in the direction of P^, and it may always be written as the centroid 

X^iT)=Xli+T-— (1) 



where M = ^—P^^P^ is the invariant mass of the system, r is the CM proper-time, both Lorentz scalars, and X^ is 
a constant 4-vector, identified as the spatial CM coordinate. Appropriately fixing the zero of t, A^ may be assumed 
orthogonal to P^ without loss of generality, Xg -P = 0. The properties of r as an observable were discussed in Refs. [26| 
and [23; those of Xq are discussed in the present article. 

With Lorentz-Poincare symmetry, the common view maintains that the CM coordinate and the internal symmetry 
of relativistic systems should be deduced from P^ and J^'' alone, together with the particles' spins, if the latter 
existii^^— . For manifestly Lorentz-covariant (observer independent) expressions for the CM coordinate, the colloquial 



choice so far has always been the so called centre- of inertia 






^ _ _:^^ (2) 

This choice was backed up by the fact that X^ is the unique solution possible for X^ if it is required to be formed 
of P'^ and J'^'^ alone^. It ceases, however, to be the solution when dependence on internal observables is taken into 
account, as is manifestly discussed in the following. 

Once the CM coordinate (Eq. ([J)) is assumed to be known, the total angular momentum may always be split into 
combination of orbital (CM) and internal parts, 

j^'^ is the internal angular momentum, relative to the centre-of-mass. From Eq. ([3]) it follows that it is a constant. 
Out of the 6 components of j'"^ 3 are independent of the CM-coordinate, fully determined by J'"^ and P^ via the 
condition 

e^uXpiJ'^'^-JnP^^O (4) 

which follows from Eq. ([3]). These components oi j^"^ constitute the spatial internal angular momentum tensor 



£'"' = A^A^ J^" = - {J'"'P^ + J^f'P'' + r^P") -^ = A^A^-^^ (5) 



where 



pfipi^ 
A"^ ^ gf^-^ + -j^ (6) 

is the 3-D metric tensor in the spatial part of the CM reference frame. £'^'^ is certainly non-zero in the general case. 
The remaining 3 components oi j^^ determine Xj^ . It is convenient to define the vector 

^ - M2 ^'> 

which incorporates these components (Q • P = 0, so Q^ contains only 3 degrees of freedom). Inverting Eq. ([3]), X^ is 
uniquely defined in terms of P^, J'^" and Q^ : 

x^o=- jj2 " =x,' + Q' (8) 

In the following, we refer as internal to observables that : 1) are invariant under uniform translations ; 2) if not 
scalars, all their components are confined to the spatial part of the CM reference frame (the hyperplane perpendicular 
to P'^). Thus, all internal vectors A^ satisfy A- P = 0, with a corresponding relation for internal tensors. 

Q^ and f^" are then constant internal quantities. By Eq. (|8]), Q^ takes the role of a shift or displacement vector, 
telling us by how much the spatial part of the CM coordinate is removed from the centre-of-inertia. It is noted that 
Q^ cannot be formed out of the Lorentz-Poincare global generators P^ and J^'^ alone, because the only internal 
quantity that they can form is i'^'^. A non-zero Q^ thus challenges the common view that the relativistic CM should 
be constructed of P'^ and J'^'^ alone^i^. Q'^, if non-zero, is therefore a new constant internal vector. 

Expressed in terms of £^'^ and Q^, j^^^ becomes 

f = t" - Q^'P" + Q'^P" (9) 

While £^'^ is responsible for (spatial) rotations relative to the CM reference frame, the remaining part of j^'^ may be 
regarded as the internal moment corresponding to the Lorentz boost. In the following we demonstrate that Q^ does 
indeed exist and is non-zero in the general case, being proportional to the LRL vector of the system, and is therefore 
responsible (together with £^^'^) for internal symmetries. Combined together via Eq. ^, j'^" is then responsible for 
the fundamental internal symmetry of the system. 



III. RELATIVISTIC INTERNAL (LAPLACE-RUNGE-LENZ) SYMMETRY 

In the present section we discuss the internal symmetry induced by the global Lorentz-Poincare symmetry. It 
has already been shown in the pasli^"— i^ that LRL vectors exist in general rotationally symmetric Newtonian-like 
systems, which generate together with internal rotations the internal symmetry. The picture will now be completed 
by showing that the mere existence of the global Lorentz-Poincare symmetry implies the existence of LRL vectors, 
which, in the same way, generate the internal symmetry together with internal rotations. 

The basis for the generalization of the LRL symmetry to general rotationally symmetric system has been showrt^^ 
to be incorporated in two basic propositions. Due to their importance, let us rephrase and prove these propositions 
in manifestly Lorentz-covariant terms. 

Even in the absence of clear and unique definition of canonical phase-space variables in classical relativistic systems, 
still the behaviour of various observables under global transformations of Minkowski space-time determines certain 
relations that any plausible definition of Poisson brackets (PB) must satisfy. These alone suffice to determine the 
generic properties of the internal Lie-Poisson algebra. Let {A, B} denote the PB of any two observables A and B. If 
6G is the generator of an infinitesimal space-time transformation, and 6A is the variation of an observable A under 
that transformation, then the PB should satisfy the relation 

5A^{A,Sg}, (10) 

together with the standard rules of Lie-Poisson algebras^, namely : 

Antisymmetry : {A,B} = -{B,A} (11a) 

Jacoby identity : {A,{B,C}} + {B,{C,A}} + {C,{A,B}} = (lib) 

Product ("Leibnitz") rule : {A, BC} = {A, B}C + {A,C}B (He) 

In the following the fundamental PB are those of the Lorentz-Poincare Lie-Poisson algebra, 

{J"", J^P} = g^'^rP - g'^^J'^P - g^Pr^ + g'^PJP^ (12) 

Besides these relations, no specific canonical structure is assumed. The PB of any other quantity which is constructed 
from PP and JP'^ are easily computed using the product rule (jllcp and the derivative rule 



{AJ{B)}^{A,B}f'{B) (13) 

which follows from it. The PB of other observables are deduced from their transformation properties. In particular, 
all 4- vectors VP satisfy 

{yP, J"^} = gP^V - gP^V^ (14) 

and all internal observables A satisfy {A,PP} = 0. 

The spatial internal angular momentum tensor iP'^ is defined from PP and JP^ via Eq. ([S]). Dual to it is the vector 
(which is proportional to the well-known Pauli-Lubanski vector) 

where IJp — PP/M is the unit 4- velocity vector of the CM-frame, with the inverse duality relation 

iP"" = eP^^PexUp (16) 

£P is also an internal quantity. The self PB of £P'^ are 

{iP'',i^P} = APPi^'' - A'^^iPP + A^PiP^ - AP^iP" (17a) 

with the corresponding ones of £P 

{ePJ'^j^eP'' (17b) 

From the rotational PB ^^, it follows that the PB of any internal vector Ap with (P'^ are 

{AP, r^} = AP^A" - AP^A^ (18a) 

for £ in tensor form, or 

{IP, A"} = {AP, t} = eP^^PAxUp (18b) 



for i in vector form, indicative of the fact that i^'^ is the generator of internal rotations in the CM reference frame. 
It is clear that the PB ([T7|) and p^ maintain the property of being internal. 

Let K^^ be an internal vector observable. K^ = K^K^ and K-£ — K'^lfj^, being internal scalar observables, certainly 
satisfy {£'', K'^} = and {£'', K ■ £} ^ 0. Then it may be shown that : 

proposition 1 The PB {K^^, K'^} are proportional to £^^" if and only if {K^^, £ ■ K} = 0. 

Proof Rf^, being an internal vector, satisfies K P — and {K^, P"} — 0. Then, by the basic PB rules (fTTj) it follows 
that the self PB {K'^,K'^} constitute an internal tensor. Being anti-symmetric, these PB may always be written as 
{Rf^, R"} = e'^'^^PAxUp with A^ some internal vector. Therefore, applying Eq. (|18bp for R^^, it follows that 

{RfJ ■ R} ^ {if^, r } R„ + {R", R""} C == {R^, R"} £u = e^"'^P£^KxUp (19) 

Since both £^ and A'' are internal vectors, it follows from Eq. ([T^ that {R^, £ ■ R} vanishes iff ^^ and A'^ are parallel. 
Thus, using the duality relation ([T6|. follows the proposition. QED 



In 2-body systems any constant internal scalar observable must be a function only of the total mass M and of 
£'^ = £'^£p = £'^'^ifj,„/2. Let R'^ be an internal vector with constant squared magnitude R^. Then necessarily R^ must 
be a function of M and £'^, say R^ = F (M, f^). For more general systems, constant internal scalar observables need 
not be functionally dependent on M and £'^ only. Still we may consider those vectors R^^ for which the product £ ■ R 
is if -invariant, {R'^, £ ■ R} = 0, and their magnitude R^ is some function R^ — F (M, £^, A) where A stands for any 
internal scalar observable which is ii'-invariant, in the sense that {R'^jA} = 0. Then we have 
proposition 2 Let R^ be an internal vector such that : 

1. The product £ ■ R is R -invariant 

2. R^ = F{M,£^,A) 

Then the self PB of Rf" satisfy 

{R^^R'^}^ ^-Q^t^" (20) 

Proof From Eq. p8b| it also follows that 

{if ^,£2} = 2{if^,r}£^ ^ 2e'"'^P£^RxUp ^ -~2£'"'R^ 
Let {is:^ , iC" } = a • £'"' . Then 

{Rf", R^] = 2 {R^', R"} R^ = 2a- £^"'R^ 
Combining the last two relations thus yields 

{R'',R^}+a-{R'',£^} =0 
from which follows, using the derivative rule (fT3|) . Eq. (|20)) . QED 

Any vector that satisfies the conditions of these propositions may be regarded as a relativistic LRL vector. In 
2-body systems there are 6 internal degrees of freedom. 4 of them are contained in the total relativistic mass M and 
the internal angular momentum £^'^ . The other 2 must be contained in a LRL vector, because it is always possible to 
construct a constant vector in the plane of motion^^ (the plane defined by £^^'^), and any constant scalar can only be 
a function of M and £'^ . 

It follows, therefore, that the shift vector Q^ must - necessarily - be a LRL vector which generates, together with 
£^^'^ , the internal symmetry of the system. This internal symmetry is governed by the PB p7)) . (|18p and ([20)1 . With 
j'^'^ being expressed in terms of £^^^ and Q^ via Eq. ^ it follows that j'"^ is the responsible for the fundamental 
internal symmetry of the system. 

As is discussed in detail in Ref. l22l the nature of the symmetry is determined by Eq. (1201) . For any particular 
system the value of 

'^--ign ^ (21) 



is the same for all LRL vectors, depending on the energetic state of the system : ij = +1 or —1 for bound or unbound 
systems, respectively. The internal symmetry generated by j'^" is then, respectively, SO{4) or 50(3, 1), for bound or 
unbound systems. The transformations generated by the LRL vector change, for a given value of the total energy, the 
internal angular momentum, thus changing the internal configuration of the system - how the particles move relative 
to the centre-of-mass - taking the system from one orbit to another, with the same energy. Explicit knowledge of the 
LRL vector may be used, as it does for classical (Newtonian) systems, to provide a full solution for the configuration 
of the system. 

Finally, we have shown that Q^^ is a LRL vector for 2-body systems. The process of CM integration may be 
performed with any number of particles, and it is indeed shown in the following (section rvT| that Q^ is a LRL vector 
for free many-body systems. The same has already been shown for post-Newtonian many-body systema^^. We may 
then conjecture that Q'^ is a LRL vector for arbitrarily large, composite relativistic systems with arbitrary internal 
interactions. 

IV. RELATIVISTIC INTEGRATION OF THE CENTRE-OF-MASS 

In the present Section we bring together the results of Sees. HI] and IIIH providing an explicit procedure for the 
computation of the shift vector Q^ and pointing at its relation with the LRL vector. 

In principle, the single particle trajectories may be parameterized each by a different time-like parameter, but for 
a common evolution picture a common parameter is required. Let a be such a common evolution parameter. Then 
the single-particle trajectories are x'^ — x'^{a). Derivatives relative to a are denoted in the following by an overdot, 
so that the particles' generalized velocities are i'^ = dx'^/da. 

Since X^' = X'^A'-^ is constant, its determining equation may be put in the form 

^-A^ = (22) 

The Lorentz-covariant generalization of the Newtonian CM, the 4-vector 

(with Mo — X)a '^a as the Newtonian total mass) does not satisfy Eq. ()22[) : however, its corresponding derivative, 

da Mo 

{v^ = i'^A(^ is the a-th particle's spatial velocity relative to the CM frame) vanishes in the non-relativistic limit 
(va/c — > 0). Therefore, the time-varying part of X^A(^ is purely relativistic. In fact, in the non-relativistic limit 
X^ A(^ — X^ — > 0. Moreover, whatever the vector X^^ may be, its behaviour under uniform translations x^ ^ x'^ + a^ 
must always be 

X^^^Xl^ + a^ + ^P^, (25) 

exactly like that of X^A^^, so the difference X^A^^^ — X^ is an internal vector. Therefore, in order to determine X^, 
we look for an internal 4-vector i?^ which satisfies 

and vanishes in the non-relativistic limit, so that X^ is given by 

Xli = X^A(; - i?'^ (27) 

The centre-of- inertia ([U, being a solution of Eq. (P^ . provides, via Eq. (I?f|) . an immediate solution of Eq. (|26l) in 
the form 

i?^ = X^A^ - Xf (28) 

If Ri was the only possible solution to Eq. (f26| . then X^ must be equal to X!^ and Q^ must vanish. This, however, 
is not the case. Non-trivial solutions, independent of R'^, are possible, with corresponding non-zero shift vector Q^ : 

Qt'^R^I;- Rt" (29) 



The process of solving Eq. ((26|) for the non-trivial solution, identifying the shift vector ((29|) and constructing 
consequently the spatial CM component X^^ via Eq. (jS]) , is referred to in the following as integration of the relativistic 
centre- of -mass. The main property of the process, namely that the time-varying part of the Newtonian CM for 
relativistic systems is purely relativistic, vanishing in the limit c — ^ oo, was first used by Dahl2^ to compute Xn in the 
post-Newtonian approximation of a 2-particle system interacting electromagnetically or gravitationally. He showed, 
via integration of an equation like Eq. ((26)) . that X^ results in a time-varying vector which is of order 1/c^ plus a 
constant of integration X^- The surprising result was that the constant X^ was not equal to Xj, the centre-of- inertia, 
but rather the difference Xo — Xi (now identified as the shift vector Q) was found to be proportional to the LRL 
vector of the corresponding Newtonian system. 

As a brief illustration of this procedure, let us recall, in an adapted way, Dahl's computation. Consider a 2- 
particle system with masses mi,m2, possible electrical charges 61,62, described in the CM reference frame with 
post-Newtonian EM/gravitational interaction^. With r the relative coordinate, v = r the relative velocity and fj, the 
Newtonian reduced mass, then in the post-Newtonian approximation Eq. (I26p becomes 



dR (mi — 777,2) 
'dt " 2M2c2 



(a'«' + 7)^ 



k{v ■ f) 



(30) 



(k — 6162 or K ~ — G77717T72 for thc electrical or gravitational case, respectively; in the CM frame R° = Q hence only 
the spatial part is relevant). Using the Newtonian equations of motion (which are sufficient since thc required overall 
accuracy is 0(1/6^)) 



dv nf 



(31) 



the square brackets in the rhs of Eq. (j30p may be expressed as a total time derivative in either of two ways 



(-^ + 7) 



k(v ■ r) 



dt 



[/i [v ■ r) v\ 



d 
'dt 



(-^ + 7) 



(32) 



The post-Newtonian trivial solution, associated with the centre-of- inertia via Eq. ([28| , is 



i?i — X^ 



X, 



7771 ~rn2 / 2 



(m-^ + 7) 



(33) 



It corresponds to the total derivative in rhs of Eq. ([32|). The other integral of Eqs. (|30l) and (|32|) combined identifies 
the non-trivial solution 






(34) 



The CM-displacement vector Q (P^j) 



Q = R\ — R2 — 



77ll — 7772 



jjLV H \r — iJL{v-r)v 

r ' 



2IvQc^ 
is clearly recognized as being proportional to the LRL vector of the corresponding Newtonian system. 



K 



(-^ + 7) 



r~-^[v-r)v^v'><.i-\ r , 

r / r 



(35) 



(36) 



(. = jifx V being the internal angular momentum vector. This is, in essence, Dahl's result^'^. 

The existence of two independent solutions to Eq. (|26p was repeatedly verified, in an analogous manner, for a pair of 
fully relativistic non-interacting particles^^; in the post- Newtonian approximation of electromagnetic or gravitational 
many-body systems-^; and for the post- Newtonian extensions of general centrally symmetric 2-body systems^. Dahl's 
results and the cited computations were all performed in the CM reference frame, but this is only a matter of 
convenience and simplicity - it may be explicitly shown that these results are valid in any reference frame, relative 
to any time-like evolution parameter, as is indeed the case in the two fully relativistic examples discussed in the 
following. And in all these systems the difference between the two independent solutions, which defines the boost's 
internal moment, is proportional to a LRL vector with a proportionality coefficient of the order of l/c^. It follows, 
therefore, that this is a generic property in relativistic systems. 



V. CM INTEGRATION FOR A TWO-BODY SYSTEM WITH SPECIAL SCALAR-VECTOR INTERACTION 

In the present Section we demonstrate the CM integration for a fuUy relativistic system with interaction. For 
simphcity and clarity of the exposition, a relatively simple system, a 2-body system with light-like antisymmetric 
interaction, is chosen : Events coupled by the interaction satisfy the light-cone condition (xi — X2) = 0, so that the 
interaction is retarded for one of the particles and advanced for the other. In this way, any event on one particle's 
trajectory is coupled, via the interaction, to a unique event on the other particle's trajectory, and a canonical structure 
is possible. Also, it is assumed that the interaction is a combination of vector (EM-like) and scalar interactions whose 
coupling constants are equal up to a sign. Defining the common coordinates 

^M — T,A' _ ^M 

such a system was discussed by Duviryakii who showed that, with the total momentum P^ canonically conjugate to 
z^ and g'' canonically conjugate to a;^, the dynamics of the system is determined by the first-class constraint 

0(.,p,.,.)^.^-2(^:^(i::^_MM!)_,(M^).o (38) 

F ■ X F ■ X 

Here 

M* - 2M2 {ml + ml) + {ml - mlf 



{AP-Ml)' 



M^ - (toi -ma)^ 



4M2 
and 



(39) 



X't 



g {M^) = ^ M^ ^ {mi - am2) (40) 

where M^ = rrii + 7712, X — sign(i:i • x) = sign (±2 • x), k, the vector coupling constant and k! = an = ±k the 
scalar coupling constant. Assuming a general evolution parameter ct, there exists some coefficient A (not necessarily 
constant) so that for any observable A(z,P, a;,^) its cr-evolution equation is 

j-A{z,P,x,q) = ^{A,(^} (41) 

with {., .} being the canonical Poisson brackets over the 16D phase-space {(z, P, x, g)}. 
The internal momentum of the system is 

n" = g^ - ^—^x^" (42) 

P • X 

This is indeed an internal vector, as verified by 11 • P = and {11^, P"} = 0. Its evolution equation is 

^ = -AAf^^< (43) 

da {P-xf 

With n^ the constraint (J3H1) is simplified to 

= n^ - -^ '- - b (Af 2) « 44) 

F ■ x 

x^^ and n'^ form an internal canonical pair, as verified by the equation 

^ = ^|^A- = An- (45) 

da 2dq'' ^ ^ 

and the fact that x^^ = {P ■ x) /M"^ . In particular, the internal angular momentum i?^" is very conveniently expressed 
in terms of x^ and 11^ as 

r" = a^a; j-^" = x'ln'' - x^n'^ (46) 



where 

is the conserved total angular momentum. Also, the centre-of-inertia is found to be 



Expressing the Newtonian CM in terms of the canonical variables, 

^ _ mix^ + 77120:2 _ 



■n^ 



(47) 



(48) 



the CM integration equation ([26l) becomes 

dm' 
da 



1 1 \ , P-x , 

T** — 2 P^ 

Af2 M2y' M4 



Af2 



Af2 



(49) 



(50) 



z^ is not an internal vector. To insure that P^ is an internal vector, we use the constancy of Xj* and obtain from 
Eq. 



P -x^ 



da \ IvP 



■W 



so that the equation for P'^ now becomes 

dRf" m- 



da 2 

An immediate solution is, of course, 



VAf2 M2 



p-x 



da V Af2 



n/^ 



7\f2 7\f2 



P-X 

1^ 



nM 



(51) 



(52) 



which is easily recognized as the trivial solution P^ = X'^A',^ — X^. To obtain the non-trivial solution it is convenient 
to transform Eq. (ISTI) . using the constraint equation (l44l) . into 



dP^ ml - m\ 



d /P-x, 



(A//2 - M^) x^-^l ^TTT^n^ 



da 2A//2Af2 ^ °^ -L dcr V Af2 

2 (7771 -7772) 6 (Af 2) 






Afo 



A/2 — (r77i — 7772)^ 
2 (7711 - 7772) 



^r^n^ 



Af„ 



Af2 — (^iTli — 1712)" 



-L dcr V Af 2 

2g(M2) 
p-x 



tfi^ 



d fP-x „ 

— n^ 

da \ A/2 



(53) 



With the help of equations (|43l) and (|45p and the relation {P-x) = APx'^, the expression within the square brackets 
is then converted into a total derivative. 



tfi^ 



2g(A/2) „ ^ ^ 2g(A/2) 



P-x 



d_ 
da 



(n-a:_L)n^- 



P-x ^ 
9 (M2) 



P-x 



Thus we are able to identify the second, non-trivial, solution of Eq. (|5ip as 



i?^ 



2 (t^i — m,2) 



A/o A/2 - (r77i - 7772)^ 



(n-a;_L)n'^ 



P-x 



P-x 



W 



(54) 



(55) 
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To compute the CM-displacement vector Q'^, it is convenient to transform the trivial solution Eq. ((52)) . using the 
constraint equation (j44l) , in a way similar to Eq. (|53p , 



< = 



(M^ - M^) a;'[ 



2M^AP 
2 (toi - ma) 6 (A'/^) 



P -x 







P-X^ 



■w- 



(56) 



Then we obtain 



i?5' - i?^ 



2 (?7ii — 7712) 



M„ 



M"^ — {mi — 7712)" 
2 (mi - m2) 



n^x^- {n-x^)w 



P-x ' 



Mo M2 - (mi - m2)^ 



,/^M 



(57) 



where 



K^" = U.I'"' 



p-x 



2^^ 



Wx' 



(n-a;_L)n^' 



p-x -L 



(58) 



is the LRL vector of the systeroii. 

The self Poisson brackets of K'^ may be computed directly using the canonicity of x^ and 11'', but this effort may 
be saved with the help of the LRL symmetry property ((20)) . using 



K^ ^b (Af 2) e^ + 



M2 



Then 



{K^',K''} 






£''•' ^-b (M^) i^"' 



(59) 



(60) 



The sign of 6(Af), which is the same as sign(Af — Mg), determines the boundness index 77 (Eq. ([11])) and thus the 
type of the symmetry. It is straight-forward to check that the PB in post-Newtonian systems (Eq.(76) of Ref. |23 ) 
is the corresponding limit of Eq. (pO]) . A remarkable feature of Eq. ([SO)) is its independence on any detail of the 
interaction, suggesting that it is universal, valid for all 2-body systems. 

Finally, we notice that since x'^ is the spatial interparticle vector r in the CM frame, the LRL vector (|58p is of 
the same structure as the Newtonian LRL vector (|36p . The similarity implies, in particular, that the orbits are fixed 
conic sections. This simplicity is due to the equality (up to sign) of the coupling constants of the scalar and vector 
interactions. To elucidate this aspect, a simplified version of this system - a particle in a Coulomb field modified by 
a scalar field - is discussed in Appendix |X| 



VI. CENTRE-OF-MASS INTEGRATION FOR NON-INTERACTING FULLY RELATIVISTIC MANY-BODY SYSTEM 



A main argument that follows from the discussion so far in the present work is that LRL vectors are derived from 
the Lorentz boost. Although the LRL symmetry is known so far to be found only in 2-body systems, the fact that any 
composite system is endowed with a Lorentz boost strongly suggests that such systems could also be endowed with 
LRL symmetry. This possibility has already been explored in the post-Newtonian approximation^^. In the present 
Section we compute the internal moment of the Lorentz boost (or, what's equivalent, the shift vector Q^) for a fully 
relativistic, non-interacting many-body system and demonstrate that it is indeed a LRL vector. This may serve as a 
starting point for establishing the relativistic LRL symmetry in general many-body systems in the future. 
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Consider a system of free, non- interacting N point particles with masses ttIq, moving on straight-line trajectories 
^ai'^a) with constant unit velocities u^ and r^ being the proper time of the a-th particle. The total linear and angular 
momenta are 



p' = J2p' 



J^^=E(^aPa-^X) 



(61) 



with the single particles' linear momenta p^ = niaU^. A generalized Lorentz factor is defined by 

/ CLTa \ 

la 



dTa 

da 



1 



V^a 



allowing us to write the particles' generalized velocities as 



i" = la'< 



Let Ea = ~Pa ■ U be the single particle energy in the CM-frame. It is convenient to introduce the notations 

{Xa ■ P) 






^ = A'^-r^ — tA" 






]VP 



IPM 



(62) 



(63) 



(64) 



for the single particle spatial coordinate in the CM reference frame, and 

<Z^^A>:^-p^-i?,t/'' (65) 

for the spatial component, in the CM reference frame, of the single particle momentum, satisfying the CM-constraint 

E'^a=0 (66) 



The single particle energies in the CM frame then become Ea — \/'(n\ + *& ^^-d the internal angular momentum is 
Substituting P^ and J^^ from Eq. (ICTI) . the centre-of-inertia ^ is 



Xf 



M2 



E 



^a ^ 



{Xg ■ P) 
Af2 



E 



T,/r ?a + 



M 



M2 



-9a 



with the trivial solution to Eq. ([26l) 

pp _ Eg "^g^^' _ yP _ V- /^!Z!2i _ :^^ fA' _ \^ (^g ■ P) „M 
1" Mo I"^Im„ Afj^" ^ M2 '^'^ 

a ^ ' a 

With the relation ([63j) . the fundamental equation of the relativistic CM integration, Eq. (pB)). becomes 



dcr 



Mo 



M„ 



(67) 



(68) 



(69) 



Using the fact that the lorentz factor also satisfies 7"^ = —Ua ■ ia, and since the particles' unit velocities {u^} and 
all the qj^ are constant in the absence of interactions, it can be shown (see Appendix [B] for details) that a non-trivial 
solution is possible only if there exists a vector G^, composed of the particles' unit velocities, that satisfies 



ETa"'«-^ = E(G •*-)'?« 



(70) 



so that 



^2 = 7^E(^-^'^)< 



Mo 



(71) 



Since Eq. ([70]) is linear in the particles' velocities i^, and these velocities are functionally independent because the 
particles' trajectories are independent, it is satisfied, again applying Eq. (p5)) . only if 



G-Ua = l 



Va 



(72) 
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Writing G^ as the linear combination G^ — 'Ylia^a'^a with constant coefficients tta, these coefficients are uniquely 
determined by the linear system 

y^ [ua ■ Ub) ab = 1 Va 

6 

thus verifying that G^ is uniquely determined by Eq. (|72p . 

The non-trivial solution (j7ip contains both spatial and time-like parts of the particles' coordinates x'^ in one term. 
To separate them, in analogy with i?^' in Eq. (|68p. we notice first that multiplying the a-th equation in ()72|) by m^, 
and summing them all yields 

G-P:=^Mo (73) 

Separating from C its CM-spatial part, 

G ■ P AT 

G^" - G'^ - P^ - G^ - i!i^pA' (7A\ 

the non-trivial solution ([71]) becomes 

^2 = ^ E (G^ • ^'^) I'a - E ^'^^^ (75) 

a a 

The trivial solution R'^ (I55|) may also be expressed in terms of the vector G^. The coefficients of ^^ may be 
expressed in terms of G^^ via the relations 

Gl_ ■ Qa _ G± -Pa _ G -Pa P ■ Pg _ TUa _ Eg , . 

Mo ^ Mo " Mo M^ ^ Mo M ^ ^ 

so that R^l becomes 

^1 = ]f^ E (G^ • la) e - E ^< (77) 

a a 

Combining Eqs. ([77]) and ([75]) for the computation of the shift vector Q^, the terms containing the time-like part of 
the particles' coordinates cancel each other and we obtain 

Q^ = ^5" - i?2^ ^ ^ E [(^^ ■ '?«) ^a - {G± ■ Ca) qH] - 
a 

1 fi my 

= Wo^^-- ^ ^'^"'^^ " ^^^'^'^ = ~ii7 (^^^ 

a 

Q^ is indeed a generalized LRL vector. For a 2-body system let q^^ = g^ = —q2- Then G^ is easily found as 

G^ = ^iH^l^H^q. (79) 

M^ — (toi — TO2) 

so that Q'' is proportional to q^f-^'^ , a more familiar form of the interaction- free 2-body LRL vector. For a many-body 
system the vector ([75]) satisfies Q^€^ = 0, and its self PB are found, taking into account that {G^ , G"} = 0, by direct 
computation, 

{Q^ Q"} = -^ {Gxi^\ G,tP] = -^r-- (80) 

On the other hand, using the identity i^'^^^x = ^^Am _ ^m^^ it follows that Q^ = \g\1^ - (G • if] /M^. From 

{G'', G""} = it follows that {Q^^, G • ^} = 0. Thus with A^ G ■ i, Q^" satisfies the conditions of proposition 2 and 
the self PB (IM|) are verified by Eq. (PH)) . This completes the proof that Q^ is indeed (proportional to) a LRL vector, 
and the many-body system enjoys LRL symmetry. 
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VII. THE SHIFT VECTOR AND CM LOCALIZABILITY 

With j'^'^ being decomposed as in Eq. (jQ)), the vector Q^ may also be regarded as generating the internal Lorentz 
boost relative to the CM reference frame. The possibility of splitting J^'', as in Eq. Q, was appreciated, at least 
in principle, from the early days of the search for the relativistic CM^^^ with various propositions for what should 
be, in the notation of the present article, the vector Q^— i^. So far in the literature an internal moment is associated 
with the Lorentz boost only if internal spin is assumed to exist as an independent entity, in which case the moment 
depends on the spin. Here, for the first time, the association of the internal Lorentz boost with the LRL vector was 
demonstrated and established. 

The preceding analysis distinguished the vectors R^ and i?2 ^^ solutions of Eq. p6)) . Although, as integrals, the 
solutions of Eq. ([26]) are defined up to an arbitrary additive constant, these particular solutions have characteristics 
that distinguish them from arbitrary integrals : Eq. ((32|) shows explicitly that Eq. ([26]) leads to two specific solutions, 
one which is proportional to the relative coordinates and another which is proportional to the relative velocities 
or momenta. The first one was identified as the vector R^ defined in Eq. (j28p . while the other was referred to as 
the non-trivial solution Rl^. This characteristic behaviour may be verified for Ri in equations p3p . ([5^ and (|68p . 
with the leading term being proportional to r, x^ and ^^ respectively, and for R2 in equations ([M]) . ([55]). and ([75]) . 
with the leading term being proportional to v, U^ and g^ respectively. This is also the case in systems discussed 
elsewhere22i2^i^. 

These two characteristic solutions define the characteristic shift vector Q^ = _R^ — i?2 , which in turn defines the 
internal moment of the Lorentz boost relative to the CM frame. A remarkable feature of this Q'^ is its being exactly 
proportional to the LRL vector of the system, thus pointing (in 2-body systems) in the direction of closest approach 
(generalized perihelion). This has been verified for various systems, both in the present paper (Eqs. (j35p and (j57p ) 
and elsewhere^^. The proportionality of Q^ to the LRL vector K^ is a very interesting aspect of the Lorentz boost, 
because Qt" could be Q'^ = aiM,i^)Kt' + l3{M,£'^)it'^K'' with arbitrary coefficients a{M,F) and li{M,P) and still 
maintain the same internal symmetry. The reason for this particular proportionality is not clear yet. 

The existence of more than one independent solution to Eq. (j26p also adds new insight into the long-standing issue 
of the localizability of the relativistic CM coordinate. With the two special independent solutions i?^ and i?2 the CM 
coordinate may be either the inertia centroid (Eqs. ^ and ^, combined) 

■W--y^ + ^- (81) 

or removed from it by the shift vector Q'', 

Jt^U p pf2 

^rW = -^^ + Q- + ^-r (82) 

Consequently, the CM coordinate cannot, in principle, be uniquely defined in a point-like manner, with Q'^ providing 
a measure of its non-uniqueness. 

Q^ is directed towards a point of closest approach. For unbound systems (collisions) there is just one such point, 
so that Q'^, in this sense, is unique. For general bound systems there is an infinite number of directions of closest 
approach, with corresponding infinite possible directions for Q^. Thus we should consider a 'centre-of-mass' domain 
relative to the centre-of-inertia, which is linear in the case of unique direction of closest approach, or circular in the 
case of multiple directions of closest approach. 

This non- uniqueness of the relativistic CM and its dependence on Q'' may be understood in the following way2^. 
The Newtonian CM coordinate Xn depends only on the particles' coordinates and masses. Thus, any value may 
be attached to it regardless the actual configuration in which the system is, and this value may remain unchanged 
even when the configuration (namely, internal energy and/or angular momentum) changes. On the other hand, 
the relativistic centre-of-inertia depends on the system's configuration via the particles' energies, and is therefore 
sensitive to any changes in the configuration which are generated by the boost's internal moment. The uniqueness 
of the Newtonian CM may therefore be regarded as reflecting its independence on the system's configuration, with 
the opposite case for the relativistic CM. The only exception is the case of equal masses, in which the centre-of-mass 
must be midway between the two particles regardless of their configuration, and indeed, as is evident from Eqs. psp . 
([57| and ([7^ . Q'' vanishes in this case. 

In Lie-Poisson algebraic terms, in analogy with quantum mechanics, localizability is synonym with requiring the 
relativistic CM coordinate to be canonical, with vanishing self PB. Assuming the fundamental PB {X^^jP^} = g'^'^ , 
it is not difficult to show that the self PB of the CM coordinates satisfy 

A';A;;{x^x''} = {g^Q-}-l-^ (83) 
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The vanishing of the rhs, implying the self PB of the shift vector Q'^ being 

'M2 



{Q^Q"^} 



(84) 



is therefore a necessary condition for X^ being localizable. However, combining Eqs. ([57]) and ([M]). the self PB of Q^ 
there is 



{Q^.Q") 



4 (mi -m2) h(M) 



Af2 



Af 2 — (ttii — 7712)^ 



.^M^ 



(85) 



and it is easily verified that the same result (in the appropriate limit) is also obtained for post-Newtonian systems. 
Moreover, the PB (j84| imply the boundness index (Eq. (|2T|)) rj = —1 which consequently implies that the system is 
necessarily unbound. Thus, even if we tried to scale Q^^ to Q'^ = a{Ri — R2) with a some appropriately chosen 
coefficient so that Q'^ satisfies Eq. (|84)). that could be possible only for unbound systems since b (M^) oc M — Mo- 
In the general case, therefore, X^^ is not expected to be canonical, thus localizable. 

This situation is somewhat similar to the one encountered with spinning point particles, even at the classical 
(non-quantum) level. It is well-known^ii^^ that such particles, endowed with conserved linear momentum p^ and 
angular momentum J'^'^ , and identifying the relativistic spin tensor as s'^'' = J'^'^ — x^p'^ + x'^p^ and the bare mass as 
"i = \/~PtiP^ J don't move on the inertia centroid (as in Eq. (jSTjl ) but rather rotate around it, thus exhibiting a helical 
motion in Minkowski space-time. The rotation is encapsulated in the uniformly rotating shift vector q^, defined in 
complete analogy with Eq. ^ as, q^ = s^^p^jm} . The difference between the two situations is that for spinning 
point particles the spin condition s^^ Xi, = is assumed, with the effect of allowing only helical orbits and limiting the 
number of independent spin components to three only, while in the preceding analysis of composite systems there is 
no such condition on the internal angular momentum j^^ (allowed to have 6 independent components), and the CM 
coordinate moves on the centroid ((M|) . shifted from the inertia centroid by the constant vector Q^ . 



VIM. CONCLUDING REMARKS 



We have demonstrated, as the main subject of the present article, that Lorentz-Poincare symmetry implies not only 
internal rotational symmetry but also the existence of LRL symmetry in composite relativistic systems. Together, 
the internal rotational symmetry and the LRL symmetry form the fundamental internal symmetry of all relativistic 
systems. 

The association of the LRL vector with the internal Lorentz boost finds its significance in the symmetries that these 
objects generate : The (global) Lorentz boost changes the state of motion ~ the way particles move - relative to the 
(external) reference frame. In a similar way, the internal boost changes the state of motion - the internal configuration 
- relative to the centre-of-mass. Spatial rotations and the Lorentz boost form together the generalized rotations in 
Minkowski space-time. Thus, the internal boost is related to the global Lorentz boost in much the same way as 
internal rotations are related to the global rotations. It is no surprise, therefore, that the boost's internal moment is 
identified with the LRL symmetry, because classically it is the LRL vector that generates the transformations that 
change the internal configuration of the systeniS^. 

To put it differently, the LRL symmetry is the internal aspect generically associated with the (global) Lorentz 
transformations, in the same way that the internal, spatial rotations are the internal aspect of the global rotations 
(the latter being defined relative to a fixed frame of reference) ; and the rotational and LRL symmetries are attached 
together internally in the same way that global rotations and Lorentz transformations form together the generalized 
rotations in Minkowski space-time. These relations are illustrated in the following diagram^^ : 



Global space-time symmetry 



Internal symmetry 



rotations 



change of 
configuration 



global 
rotations 

X 

internal 
rotations 



-h 



Lorentz 
transformations 

X 

LRL 



The LRL symmetry is therefore geometrical in nature, as the LRL vector is associated with the internal moment 
(relative to the centre-of-mass) corresponding to the Lorentz boost. LRL vectors and LRL symmetry are therefore 
universal, characteristic of all relativistic systems. 
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In conclusion, Lorentz-Poincare symmetry demands the extended internal symmetry, manifested via the corre- 
sponding internal moments and closely associated with the determination of the relativistic centre-of-mass, providing 
us with new tools for the investigation of the internal dynamics of relativistic systems. As a consequence, the CM 
coordinate of composite relativistic systems should be explicitly constructed not only of the global quantities P^ and 
J^"^ , but also of the LRL internal dynamics of the system. 

Appendix A: The relativistic scalar-Coulomb case with infinite central mass 

It was noted at the end of section |V] that the LRL vector (|58)) is of the same structure as the Newtonian LRL vector 
(|36p . This simplicity is due to the equality (up to sign) of the coupling constants of the scalar and vector interactions. 
To elucidate this aspect, let us consider a relativistic particle with variable mass 

^\r)^m^+'^ (Al) 

in a Coulomb field with fixed centre (infinite central mass), with the Hamiltonian 



H = \lp^ +'m?+'^ + - (A2) 

in the CM frame. In a configuration with given energy H — E and internal angular momentum i ^ r x p, we isolate 
the squared momentum as 



f / K\2 , k'2 ^ ^ 2kE 



p'-pI + Z2 = [E--) -m'--^^E'-n,'-^^~'—:^ (A3) 



Equating now, up to a sign, the coupling constants, k = ±k', results in the canceling of the 1/r^ terms on the rhs of 
Eq. (jA3l) . keeping the centrifugal term £^/r^ as in Newtonian mechanics and leading to the equation 



1 drV n kEV {E^-m^)e + E^K^ 



^ ■ dO ) \r P ) £4 

with the solution 



(A4) 



1 kE J{E^~m^)P + E'^K^ 

- ^ ^ ' cos 61 (A5) 



r 



P P 



The solution is a fixed, non-rotating, conical curve, as for the non-relativistic Coulomb case. 

The similarity in form of the solutions implies similarity of the LRL vectors, which is in the present case 



It is a constant vector with magnitude 



K 



K V 

K = ex p f (A6) 

r 



= {E^ - m^) e^ + E^K^ , (A7) 



directed along the major axis of the conic section. It would be interesting to study the quantization of this vector, 
because, unlike the non-relativistic case, here it depends on energy eigenvalues, both positive and negative. 

Evidently, if k ^ ±k' then the extra l/r^ term on the rhs of Eq. (|A3|) would modify the centrifugal term and lead to 
a rotating conic section with the corresponding complicated LRL vector, as in the relativistic Kepler- Coulomb system 
with fixed centrei^i^i^S, or the Newtonian Kepler- Coulomb system modified with 1/r^ potential^!. 

Appendix B: Proof of the uniqueness of the solution Eq. I|71|l 



In the following we prove that Eq. ((70| is necessary for an internal non-trivial solution Rt^. In an A^-body system, 
let 1 < TO < A. Then, by means of the CM-constraint ([M)) . g^ is eliminated from the sum '^^Ja^la i^^ ^Iq. (IM|) to 
yield 

Y.^a'qii = Y.i"'a'-i;-n')qii (bi) 
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Since all the A*"— 1 remaining q^'s are dynamically independent, and since the sum J^a'ya ^1a must lead to an internal 
integral, then, taking into account the identity 7"^ = —Ua ■ ia, it follows that for each of the coefhcients in Eq. (|Bip 



there must be a vector Ga , independent of the particles' coordinates and depending only on the particles' unit 
velocities -u^, so that 7^^ — 7,7^^ = Ga ■ [ia — im)- 
The use of the particular index m in the resultant sum 






^m) Qa 



is of course arbitrary, and another index 1 < n < N,n ^ ni could be used instead with the sum 



E 



Eliminating g,'^j from the sum in Eq. (jBSp . the latter becomes after some algebra 

-j_r;(") . a - r )at^ 



(B2) 



(B3) 



Y.^a'€= E 



€ 



(B4) 



The sums in Eqs. (IB2p and (|B4p . must be identical, but their identity can be realized iff all the vectors Ga 

are identical, Ga = G^ \/ni, a. This completes the proof of the uniqueness of the form (|7ip for the non-trivial 
solution. 



^M. H. L. Pryce, "The mass-centre in the restricted theory of relativity and its connexion with the quantum theory of elementary 

particles", Proc. R. Soc. London, Ser. A 195, 62 (1948). 
^T. D. Newton and E. P. Wigner, "Localized states for elementary systems", Rev. Mod. Phys. 21,400-6 (1949) 
^G. N. Fleming, "Covariant position operators, spin and locality", Phys. Rev. 137B, 188—197 (1965) 
*D. Alba, L. Lusanna, and M. Pauri, "Centers of mass and rotational kinematics for the relativistic n-body problem in the rest-frame 

instant form", J. Math. Phys. 43, 1677-1727 (2002) 
^D. Alba, H. W. Crater, and L. Lusanna, "Hamiltonian relativistic two-body problem: center of mass and orbit reconstruction", J. Phys. 
A; Math. Theor. 40, 9585-607 (2007) 

H. Goldstein, C. Poole and J. Safko, Classical Mechanics (Addison- Wesley, New-york, 2000). 
'^H. V. Mcintosh, "Symmetry and degeneracy", in Group Theory and its Applications, Vol. II, edited by E. M. Loebl (Academic Press, 

New York, 1971) pp. 75-144. 
^H. Bacry, H. Ruegg and J. Souriau, "Dynamical groups and spherical potentials in classical mechanics", Coinm. Math Phys. 3, 323-333 

(1966). 
®D. M. Fradkin, "Existence of the dynamical symmetries 04 and su:j for all classical central potential problems". Prog. Theor. Phys. 37, 

798-812 (1967). 
^''N. Mukunda, "Dynamical symmetries and classical mechanics", Phys. Rev. 155, 1383-6 (1967). 
^^A. Peres, "A classical constant of motion with discontinuities", J. Phys. A: Math. Gen. 12, 1711-3 (1979). 
A. Holas and N. H. March, "A generalization of the Runge-Lenz constant of classical motion in a central potential", J. Phys. A: Math. 
Gen. 23, 735-49 (1990). 
^^T. Yoshida, "Rotating Laplace- Runge-Lenz vector leading to 2 relativistic Keplers equations", Phys. Rev. A 38, 19-25 (1988). 
Argiieso and J. L. Sanz, "Post-Newtonian extensions of the Runge-Lenz vector", J. Math. Phys. 25, 2935-2938 (1984). 
DrozVincent and P. Nurowski, "Symmetries in predictive relativistic mechanics", J. Math. Phys. 31, 2393-8 (1990). 
P. Horwitz, "Dynamical group of the relativistic Kepler problem", J. Math. Phys. 34, 645-8 (1993). 

Duviryak, "Symmetries of the relativistic two-particle model with scalar-vector interaction". Nonlinear Math. Phys. 3, 3728 (1996). 
H. Johnson and B. A. Lippmann, "Relativistic Kepler problem", Phys. Rev. 78, 329 (1950) 
C. Biedenharn, "Remarks on the relativistic Kepler problem", Phys. Rev. 126, 845-51 (1962) 

T. Khachidze and A. A. Khelashvili, "The hidden symmetry of the Coulomb problem in relativistic quantum mechanics: From Pauli 
to Dirac", Am. J. Phys. 74, 628-32 (2006) 
^-"^T. T. Khachidze and A. A. Khelashvili, Dynamical Symmetry of the Kepler- Coulomb Problem in Classical and Quantum Mechanics 
(Nova Science, New- York, 2008). 

Ben-Ya'acov, "Laplace- Runge-Lenz symmetry in general rotationally symmetric systems", J. Math. Phys. 51, 122902 (2010) 
P. Dahl, "Physical origin of the Runge-Lenz vector", J. Phys. A: Math. Gen. 30, 6831-6840 (1997). 
P. Dahl, "Physical interpretation of the Runge-Lenz vector", Phys. Lett. 27A, 62-3 (1968). 

Ben-Ya'acov, "Localization of the relativistic centre-of-mass from internal dynamics", J. Phys. A: Math. Theor. 42, 225401 (2009) 
Ben-Ya'acov, "Internal time observable of classical relativistic system", J. Physics A : Math. Gen. 39, 667-683 (2006) 
Ben-Ya'acov, "Internal-time and dilatations in classical relativity", J. Physics: Conference Series 66, 012009 (2007) 
Ben-Ya'acov, " Laplace- Runge-Lenz-like new constant in many-body systems from post-Newtonian dynamics", J. Phys. A: Math. 



14p^ 
15p 
16L. 

"A. 

18m 

19L. 
20 rp 



22u 
23j. 
24j. 

25u 
26u 

27u 
2Su 



Theor 42, 375210 (2009). 



17 

^®M. Lorentc and P. Roman, "General expressions for the position and spin operators of rclativistic systems", J. Math. Phys. 15, 70-74 

(1974) 
^"L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields (Pergamon Press, Oxford, 1975). 
^^H. C. corben. Classical and Quantum Theories of Spinning Particles (Holden-Bay, San francisco, 1968). 
^^N. Kudryashova and Y. N. Obukhov, "On the dynamics of classical particles with spin", Phys. Lett. A 374, 3801-5 (2010) 
^^J. E. Marsden and T. S. Ratiu Introduction to Mechanics and Symmetry (Springer, New- York, 1999). 
^^W. H. Hcintz, "Runge-Lenz vector for nonrelativistic Kepler motion modified by an inverse cube force". Am. J. Phys. 44, 687-94 (1976). 



